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Airframe and Control System Variables 
vertical velocity of airplane center of gravity 
vertical velocity command from autopilot 
vertical acceleration of airplane center of gravity 
gross weight 

airplane moments of Inertia about body axes 

airplane product of inertia with respect to x and z body axes 

Laplace transform variable 

airspeed 

airplane longitudinal inertial velocity 

airplane longitudinal inertial acceleration 

airplane longitudinal acceleration command from autopilot 

localizer error 

airplane angle of attack 

airplane angle of sideslip 

control system input to ailerons 

control system input to elevator 

pilot control input to vertical controller 

pilot control Inputs to attitude control system 

control system input to rudder 

control system inputs to power lever system 

power lever system input to fuel control 

power lever system input to blade control 

pilot control input to horizontal acceleration controller 

roll control system input to blade control 
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6 0 pitch control system input to blade control 

6^ yaw control system input to louver control 

6^ feedback quantity for roll controller derived from blade actuator 

position 

Sq feedback quantity for pitch controller derived from blade actuator 

position 

feedback quantity for height controller derived from blade 
actuator position 

Ah altitude error 

<(> roll angle 

$ first time derivative of $ 

$ second time derivative of <p 

<t> c roll angle commanded by autopilot 

6 pitch angle 

6 first time derivative of 0 

6 second time derivative of 6 

0 C pitch angle commanded by autopilot 

thrust vector angle (90° normal to body x axis) 

C_ T thrust vector command input to louvers 
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Propulsion System Parameters 

General form: ABC x 1 x 2 y> where ABC refers to the physical quantity, 

x x and x 2 refer to the location within the engine, and y refers to the 
engine number (see figs. 1 and 2). 

ALT altitude 

BETAF fan blade angle 

BETN command to fan blade actuator 

DB fan blade actuator deadband 

DTT ambient atmospheric temperature 
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FGROS 

HPF 

HP1T 

HP2T 

HP2C 

I HP ,I LP 

MNO 

NH 

NL 

PCN4D 

PCNHR 

PS 

PT 

SM 

TT 

V 

U 

WA 

WFH 

WG 

XNPCT 


net thrust 

gross thrust 

power absorbed by fan 

engine power available to fans 

engine power available to compressor 

power absorbed by compressor 

moments of inertia of high- and low-pressure engine rotors 
Mach number at inlet 
engine angular velocity 
fan angular velocity 

engine corrected speed command, percent 

engine corrected speed, percent 

static pressure 

total pressure 

fan stall margin 

total temperatures 

flow velocities 

engine main flows 

duct mass flows 

fuel mass flow rate 

mass flow of combined fan and engine nozzles 

fan mechanical speed, percent 

time constant in heat-soak transfer function 
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SIMULATION STUDY OF THE INTERACTION BETWEEN THE PROPULSION AND 


FLIGHT CONTROL SYSTEMS OF A SUBSONIC LIFT FAN VTOL 

Bruce E. Tinling 
Ames Research Center 

and 

Gary L. Cole 
Lewis Research Center 


SUMMARY 


The possibility of interactions between the propulsion and flight control 
systems of a three-fan subsonic VTOL aircraft has been studied using nonreal- 
time simulation. The objectives were to obtain time histories of critical 
internal engine parameters and to identify and analyze possible deleterious 
effects of engine dynamics on flight control. Analysis of the engine data was 
beyond the scope of the study. 

The simulation results indicated no deleterious effects of engine dynamics 
on the control of the aircraft with the exception of the effects of the fan 
actuator deadband. A method of alleviating these effects through feedback of 
the actuator output to the flight controller was developed. 


INTRODUCTION 


During the critical transition and hovering flight phases, the propulsion 
system of a VTOL aircraft must provide control forces and moments as well as 
the forces necessary to sustain flight. Consequently, dynamic Interactions 
between the propulsion and flight control systems are a possibility because 
logic to enhance handling qualities requires closed-loop control in which the 
propulsion system is an element of one or more high-gain control loops. 
Furthermore, the propulsion system must accommodate unique requirements such 
as power management between multiple engines and compensation for engine 
failures. 

The handling qualities of VTOL aircraft with advanced controllers have 
been studied intensively, including real-time piloted moving-base simulations 
(see, e.g., refs. 1 and 2). In these studies, the propulsion system was 
approximated by a simple input /output relationship that represented the 
response to control inputs. This permitted an evaluation of the handling 
qualities in the powered-lift mode but did not provide information on the 
possibility of exceeding Internal engine limits due to control activity. 
Further, the simple engine representation eliminates many possibilities of 
deleterious interactions between the dynamics of the engine and its control 
system and the flight control system. It is the purpose of the study to 



examine these interactions for a subsonic VTOL aircraft with lift-cruise fans. 
The objective is to provide data from which it can be determined if a high- 
gain flight control system performs satisfactorily when the Internal dynamics 
of the engine are represented and if control activity has an important Impact 
on the engine design or life. Examples of the latter are excessive turbine 
temperature cycling and fan blade actuator activity. 

The aircraft chosen for the study is a subsonic, three-fan VTOL design 
proposed as a research and technology aircraft. This design has been studied 
extensively in piloted simulations at Ames Research Center (see refs. 1 and 3). 
The propulsion system consisted of two turbofan engines, one turboshaft 
engine, and one remote lift fan. The fans were all connected with shafting 
through a gear box as illustrated in figure 1. The propulsion model was suf- 
ficiently detailed to permit monitoring of internal engine variables during 
transient operation. It Included dynamic effects such as heat soak, rotor 
inertia, core speed control, and fan blade actuation. 

The flight control system chosen for this study was the State Rate Feed- 
back Implicit Model Follower described in reference 2. This controller pro- 
vides an input /output relationship that approximates any selected second-order 
system; it provides good gust alleviation and cross-axis decoupling and is 
self trimming. This controller was chosen because it provides satisfactory 
handling qualities throughout a precision powered-lift instrument approach and 
in hovering flight. This performance is obtained through high-gain control 
loops that include the engine. Combining this controller with the detailed- 
engine model is therefore likely to expose dynamic interaction between flight 
control and engine dynamics existing for this type of aircraft. 

The study was performed in non real- time digital simulation. The aircraft 
was required to follow a specified flightpath under automatic control in the 
presence of turbulence, initial vertical and lateral offsets, and engine fail- 
ures. Example data are presented that illustrate the performance of the flight 
control system and the variation of critical engine parameters. A list of all 
of the flight conditions for the simulated approaches is Included. The time 
histories of relevant parameters for these approaches have been recorded on 
magnetic tape and are available to agencies of the U.S. Government and its 
contractors upon request. 


SIMULATION MODEL 


The Navy research and technology aircraft was chosen for this study for 
several reasons. First, it represents a recent VTOL aircraft design for which 
a complete documented set of aerodynamic and design data is available (see 
ref. 3). Second, a digital computer program used for piloted simulation at 
Ames Research Center was readily available. Finally, and perhaps most impor- 
tant, a detailed simulation of a shaft-driven lift-fan propulsion system was 
available for this aircraft. This simulation was developed at Lewis Research 
Center. 
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Airframe 


The aircraft (fig. 1 and table 1) was representative of a research air- 
craft with a wing span of 13.54 m and a gross weight of about 13,000 kg. 

Control in the powered-lift mode was provided through variable fan blade pitch, 
the thrust vectoring hoods and louvers, and attitude (see figs. 3 and 4). 

Pitch and roll control were provided through the variable pitch fans; yaw con- 
trol through differential thrust vectoring; and lateral velocity control 
through roll of the aircraft. Lateral velocity control for low-speed flight 
was also available through thrust vectoring, but this mode was not used in 
this study. 


Attitude and Flightpath Controllers 

The attitude and flightpath controllers chosen utilize the principles of 
the state rate feedback implicit mode4 following controllers described in 
reference 2. This approach was chosen for several reasons. First, it has been 
found to provide better handling qualities for the lift-fan aircraft than 
others which have been studied in simulation of VTOL aircraft at Ames Research 
Center (see ref. 2). Second, the high-gain controllers should require rapid 
response from the propulsion system, thereby exposing deleterious interactions. 
The controllers chosen for the study provided attitude command in pitch and 
roll, rate command with turn coordination in yaw (washed out at speeds below 
20 knots) , and vertical velocity and longitudinal acceleration command for 
control of the flightpath in the vertical plane. As previously noted, lateral 
flightpath control was provided through roll. 

A curved decelerating approach was specified for the study, and command 
signals necessary to follow the trajectory were provided by a flight director 
directly coupled to the controllers. The elements of this flight director are 
described in reference 2. Signals were generated to command the roll angle 
for lateral flightpath control, and the attitude rate and longitudinal accel- 
eration for longitudinal flightpath control. The pitch attitude was commanded 
to be zero degrees, which is considered to be the best touchdown attitude, for 
the entire trajectory. 

Block diagrams of the state rate feedback implicit model following con- 
trollers are shown in figures 5 and 6. These controllers differ in detail from 
those shown in reference 2 in that the manual control provisions and alterna- 
tive modes of a operation have been deleted. In addition, the controllers for 
$, 6, h, and V x have been modified to provide separate control signals to the 
aerodynamic and to the propulsive controls. Also, in those instances where 
control of the variable pitch fans is required, provision has been made for 
feedback from the actuators to the controller. The reasons for these latter 
modifications to the controllers will be developed in the discussion of the 
results. 

Details of the power lever system, the engine speed control, and the fan 
blade actuation system are shown in figures 7, 8, and 9. These systems are 
similar to those employed in the piloted simulation study reported in refer- 
ence l with the exception of the fan blade pitch actuators. These actuators 
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were assumed to have a deadband as illustrated in the lower part of figure 8. 
The magnitude of the deadband was nominally set at 0.5" and the actuator time 
constant at 0.1 sec. 

The actuators for all aerodynamic controls had no deadbands and were 
modeled as described in reference 3. The tine constant for these actuators 
was typically 0.1 sec. 

The flight director and various elements of the flight and engine control 
systems require data on attitudes, rates, accelerations, engine speeds, and 
so on. For this study, all measurements to provide these data were assumed to 
be perfect. 


Propulsion System 

As indicated in the sketch of the aircraft in figure 1, the propulsion 
system consists of two turbofan engines, one turboshaft engine and one lift 
fan. The units are all connected with shafting through a gear box. The remote 
unit is disengaged by clutch during conventional flight. The core engines are 
modified Detroit Diesel Allison XT701-AD-700 engines and the fans are Hamilton 
Standard 1.57 m (62 in.) diameter variable-pitch, low-pressure ratio fans. 
Additional information concerning these units is given in reference A. 

The propulsion representation used in this model is nonlinear and con- 
siderably more detailed than that originally used for piloted real-time simu- 
lation. The major advantage of the detailed model is that it allows monitoring 
of internal engine variables (i.e., temperatures and pressures throughout the 
engine) during transient operation. The propulsion model includes dynamic 
representations of the rotor inertias, heat soak, the core speed controller 
and the fuel valve. Block diagrams of the propulsion system and core speed 
control are given in figures 2 and 8, respectively. Each propulsion unit has 
its own representation. However, only two representations are shown in fig- 
ure 2 since units 1 and 2 are identical. The fuel controls for all three 
engines are identical. 

The major inputs to the propulsion system are fuel flow (WFM) to each 
engine (output of the speed control and fuel valve of figure 8 which will be 
discussed later) and the fan pitch angles Si, 82, and 83 (see fig. 9 ). Air- 
craft altitude ALT la input and is used to determine ambient temperature and 
pressure from standard atmosphere tables. Hot or cold days can be simulated 
by the input DTT (fig. 2 ). Ambient conditions (P 0 , TO) plus the aircraft 
Mach number MNO are used to compute free-stream total conditions from the 
compressible flow equations. The inlet pressure recovery was assumed to be 
1002 for all three inlets; therefore total conditions into the fan are the 
same as free-stream totals. Also input to the fan are rotor speed, NL, and 
fan pitch angle, BETAF. 

The fans are represented by maps of total fan corrected airflow, pressure 
ratio and temperature ratio as functions of fan pitch angle (BETAF) , and fan 
corrected speed. The maps are constructed according to the assumption that 
each fan operates along a single operating line depending on BETAF. The 
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fan maps cover the following ranges; fan corrected speed, 70 to 110%; and fan 
pitch angle, -20 to +7.3°. 

Separate temperature and pressure ratios are calculated for the hub (sta- 
tion 25) and tip (station 13) of fans 1 and 2. Inputs to the compressor are 
total pressure and temperature (PT25, TT25) from the fan and high-pressure 
rotor speed NH. The compressor is represented as single operating lines of 
efficiency and corrected airflow as functions of corrected speed. Values are 
mapped for corrected speed in the range of 65 to 107.5%. Actual compressor 
airflow W25 was calculated and subtracted from total fan air W1 (propulsion 
units 1 and 2 only) to give the fan duct airflow (WA13). 

The fan duct is simply treated as a constant total pressure loss resulting 
in PT18. Inputs to the combustor are pressure PT3, temperature TT3 (which 
is the output of compressor heat-soak lead-lag dynamics), airflow W3, fuel 
flow WFM (output of the fuel valve, fig. 8) and the compressor exit static 
pressure PS3. The combustor generates a temperature TT4P and a function of 
airflow W3, fuel flow WFM and temperature TT3. The high-pressure turbine 
inlet temperature TT4 is the output of turbine heat-soak lead-lag dynamics, 
and PT4 is assumed to be equal to PS3. The calculation of PS3 requires 
an iteration process. Pressure PT42 and temperature TT42 out of the high- 
pressure turbine are assumed to be directly proportional to PT4 and TT4, 
respectively. Temperature, pressure, and flow at the nozzle for units 1 and 2 
are calculated assuming mixing of the duct and low-pressure turbine airstream. 

The gross thrust FGROSS Is calculated from conservation of momentum 
using WG18 and VI8 (which depends on TT18 and PT18). In the case of propul- 
sion units 1 and 2, the nozzle air temperature is higher than that from the 
fan duct alone because of mixing with the hot core air. The higher temperature 
results in higher thrust than would be obtained from airflow at the lower fan 
duct air temperature. In the case of propulsion unit 3, only the fan is 
assumed to produce thrust (cure thrust is negligible). The gross thrusts are 
output to the aircraft simulation along with inlet airflows for the calculation 
of net thrusts. 

Rotor speeds NH and NL are determined from the power absorbed by the 
compressor HP2C or the fan HPF and the power output by the corresponding 
turbine (HP2T or hPIT). n the case of the fans, since all three are con- 
nected by shafting, the HPF's and HPIT's are summed as shown in figure 2. 

The time rate of change of speed NL is directly proportional to the AHP's 
and inversely proportional to rotor inertia 3I^p and rotor speed NL. NL 
is integrated to obtain NL which is the same for all three fans. Although 
not shown in figure 2, an option is available to allow declutching of the nose 
fan (no. 3). In that case, the power absorbed by fan 3 HPF3 and the inertia 
of fan 3 go to zero. 

High-pressure rotor (compressor) speed NH is determined in the same 
manner except that each rotor works independently. Therefore, NH's are not 
necessarily the same. For simulation of a core failure, the power output of 
the associated low-pressure turbine HPIT goes to zero. Also, the tempera- 
ture of the nozzle flow is set equal to the fan duct air temperature (i.e., no 
heating from core). 
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Core Speed Control Details 


A block diagram of the engine core speed control is shown in figure 8. 
Included are the dynamics of the compressor-inlet total- temperature (TT25) 
sensor and the fuel metering valve. Both are first-order lags. The inputs t • 
the fuel control are sensed compressor inlet temperature T25SN, compressor 
exit static pressure PS3 (i.e. , the sensor is assumed to be perfect — no 
dynamics), high-pressure rotor mechanical .speed NH and demanded high-pressure 
rotor corrected speed PCNHD. The demanded speed for each engine is determined 
from the fuel flow demanded by the power lever system by an equation of the 
form PCNHD =» Ci*(PCNHC)**0.2 + Cj (C i and C 2 are constants but different for 
different engines) (see upper part of fig. 8). The fuel control incorporates 
a proportional-plus-integral controller and acceleration and deceleration 
schedules. The control limits engine operation between mechanical speeds of 
57% and 1 10% . The output of the fuel valve WFM goes to the combustor as 
shown in figure 2. 


SIMULATED FLIGHT CONDITIONS 


The flight conditions for the digital simulation were specified to permit 
study of the following on the response of the aircraft and its propulsion 
system; 

* Initial localizer and attitude errors 

* Engine failures 

* Gross weight (minimum 11,000 kg; maximum 12,960 kg) 

* Ambient temperature (standard day 15°C; hot day 32.2°C) 

* Turbulence level (0.25 m/sec to 2 m/s ec rms) 

* Fan blade angle deadband 

* Fan blade actuation time constant 

■ Methods for alleviating effects of fan blade actuation deadband 
A complete list of the simulated approaches is given in table 2, 

For all of the data, the same commanded trajectory and wind conditions 
were specified. The trajectory consisted of a curved decelerating approach 
that was initiated at an altitude of about 460 m on a three-degree glideslope. 
The commanded deceleration was constant except tor the first few seconds and 
the last few seconds prior to hover. The data records taken consist of approx- 
imately the last 120 sec of the approach. The heading was 90* in all cases 
with a 20-knot mean wind from 60*. Turbulence was introduced for all caset . 
since it was found that under some conditions a small disturbance was required 
to Induce attitude oscillations caused by the fan blade actuator deadband. 

The lowest level of turbulence used was 0.25 m/sec rms. The wind model used 
and the Dryden turbulence model are defined in MIL-F-8785B (see ref. 5). The 
random number sequence used in generating the turbulence was identical for 
all trajectories so that the results sre directly comparable. 
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RESULTS AND DISCUSSION 


Typical Results 

An example of typical results obtained from the simulation is given in 
figure 10. For this set of data, the gross weight was set at the maximum 
(12,960 kg), the ambient temperature corresponded to a hot day (32.2°C), a 
20-knot wind was imposed at an angle of 30° from the flightpath with an rms 
turbulence level of 1 m/sec, and an engine failure was introduced at an atti- 
tude of about 305 m. In addition, the initial altitude and lateral positions 
were displaced 20 ra from the commanded flightpath. Fcr these rest'ts the 
actuator deadband was zero. The results indicate excellent flightpath track- 
ing with an altitude loss at the time of the engine failure of only 0.6 m. 

The time histories of the internal propulsion system parameters show the 
dynamic response to turbulence and to engine failure. Not; that, with one 
engine inoperable, the peak temperature at the burner outlet of the operating 
engines reaches about 1620 K when velocities approaching hover are reached. 

This is the temperature specified for the one-hour contingency rating in the 
engine design (see ref. 6). 

A complete list of conditions for which similar data were obtained is 
given in table 2. These results do not include the final deceleration to 
hover. Results that are limited to approximately the last 20 sec prior to 
hover and which include the effects of temperature rise due to reingestion are 
discussed in a subsequent section. A complete list of the conditions for which 
these results were obtained is given in table 3. 


Effects of Fan Blade Actuator Deadband 

A characteristic of hydraulic actuators is that some deadband will exist* 
the magnitude depending on the particular design. A block diagram represent- 
ing the dynamics of the fan blade pitch actuator with its deadband is shown in 
the lower part of figure 9. Note that the output nf the actuator is an input 
to the engine model, and therefore the actuator deadband can be easily included 
in simulations using either a detailed engine model or a simpler representa- 
tion such as used for piloted simulation (ref. 3). The effective actuation 
time constant caused by the deadband will depend on the magnitude of the com- 
mand signal, being relatively small for large Inputs, and becoming infinite 
for Infinitesimal signal inputs. For this study, a nominal deadband of 1/2* 
was assumed. 

The effect of the deadband on the angular acceleration and on the fan 
blade pitch actuation is shown in figure 11. For this approach, mild turbu- 
lence was imposed to provide a continuous disturbance. It would be antici- 
pated from the nature of the turbulence model and from the decreasing flight 
speed that the envelope of the angular accelerations would decrease as the 
approach continued, as was the case with zero deadband shown in the upper 
right part of figure 10b. When the fan blade actuation deadband is present 
and the airspeed has decreased to about 60 knots (tims *70 sec), a high- 
frequency oscillation in attitude occurs that is caused by a limit-cycle 
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oscillation of the fan blade pitch actuators. This activity can be more 
clearly seen in figure lib, where the actuatoL command signal and actuator 
position are shown on an expanded scale. This oscillation does not cause any 
appreciable deviation in the commanded attitude or flightpath but would cer- 
tainly be annoying to a pilot, might induce resonant oscillations if the air- 
craft structure, and would decrease the service life of the actuators. 

The limit-cycle oscillation is a consequence of the low airspeed and the 
attendant decrease in aerodynamic damping. This damping must be provided by 
the control system, requiring the fan blade actuators to follow small inputs. 
This limit-cycle oscillation can be demonstrated for a much simpler system 
than the lift-fan aircraft and its control 'vstem (see fig. 12). The con- 
trolled element in this case is a first-cru^r rate damped system representa- 
tive, for example, of single degree of freedom roll control of a conventional 
airplane. The controller is a simple state rate feedback implicit model fol- 
lower developed using the same principles described in reference 2. For suf- 
ficiently high values of the gain K, the overall transfer function will be 
that of a simple rate damped system with a time constant 1 /K i regardless of 
the value of the damping parameter Di, providing that the actuator dynamics 
can be approximated as l/(t e S +1). However, when a deadband is present in 
the actuator of the type assumed for the lift-fan aircraft and the damping of 
the plant is low, the response of the system to a small pulse consists of a 
limit-cycle oscillation as illustrated in figure 12. 

It was shown in reference 2 that the path 1-2 in figure 12 must contain ^ 
dynamic element to be physically realizable. Further, it was shown for a 
linear system that representation of the actuator dynamics in the path 1-2 
yields the desired system response. This result leads heuristically to an 
attempt to regain the desired performance for a system with a deadband in the 
actuator by placing a model of the actuator, inducing the deadband in the 
path 1-2. The results of such an experiment are shown in figure 13. It is 
shown that this approach eliminates the limit-cycle oscillation, although the 
response differs from the desired linear transfer function represented by 
x/x c ^ Ki/(s + Ki), However, for inputs equal or larger in magnitude to the 
deadband little difference from the desired linear response can be observed. 


Propuls ion /Control System Integration to Suppress the Effects 
of Fan Blade Pitch Actuator Deadbands 

Application of the heuristic approach to suppress the effects of the 
actuator deadbands is more complicated than for the simple system described in 
the preceding section. The controller , as developed in reference 2, controls 
both the aerodynamic and the propulsive forces and moments* Application of 
the technique Indicated by the study of the simple system requires that a 
model of the actuator be In the controller for each force or moment generator; 
otherwise , limit-cycle tendencies were found to exist because the activity of 
all actuators was not properly modelled. Th* requirement leads to developing 
the roll, pitch, and vertical flightpath controllers : nto branches (figs. 5 
and 6) so that appropriate actuator dynamics can be insisted in each. Further, 
there is a possibility that the commands to the aerodynamic and propulsive 
actuators will diverge, yielding opposing control moments. This is prevented 
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in the long term through an interconnect between the branches which has a suf- 
ficiently long time constant so as not to interfere with the short period 
dynamics. When the airspeed becomes sufficiently low, it is necessary to 
eliminate inputs to the aerodynamic controls since their activity further 
excites limit-cycle oscillations at these speeds. 

The success of the above modifications in eliminating the limit-cycle 
oscillations can be described as erratic. There is nothing in the mechaniza- 
tion of the actuator model to guarantee that it is operating synchronously 
with the actuators. Accordingly, under some conditions, the mechanization 
reduced the level of the oscillation; under others, the oscillations increased. 
One of the latter cases is shown in figure 14. It can be seen that, compared 
to the uncompensated case, the pitch frequency is reduced and the amplitude is 
increased over the final 20 sec of the approach. These two changes are likely 
to be deleterious to piloted operation since there is sufficient time during 
an oscillation to reach attitude errors of the order of a degree, and the fre- 
quency is sufficiently low for the pilot to attempt corrective control — 
possibly leading to pilot-induced oscillations. 

For the simple single degree of freedom system discussed earlier, the 
obvious solution in overcoming the deficiencies of the actuator model approach 
is to remove the actuator model from the controller and replace it with the 
output signal from the actuator itself multiplied by 1/K (see fig. 12). This 
ensures that the input to the controller and the actuator are synchronized. 

An additional advantage is that the actuator is likely to require a far more 
complex representation than the simp. a form shown in figure 9, and accurate 
modeling might be difficult or at least an unwanted complication. The feed- 
back of the actuator output to the controller eliminates this complication as 
well as ensures that the controller is receiving correct information. 

For the aircraft controller, the solution is not as straightforward. As 
indicated in figure 9, the signals are mixed, each actuator being driven by a 
linear combination of the signals from each controller, that is: 


BETN1 


S 

BETN2 

- [M] 

6 0 

BETN3 




when H is the matrix of gains in the signal mixing logic. Estimates of the 
quantities to be fed back that indicate the effects of actuator lags and dead- 
bands were obtained by multiplying the inverse of M by the output of the 
actuators; that is. 


K ' , S* t e 


BETAF1 

K e«9 

- [M -1 ] 

BETAF2 

©*> 

-o* 


BETAF3 
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The quantities K v ^, Kq, and represent the gains between the input to the 

signal mixer and the points at which 6 t g, 6 0 , and 6^ are fed back to the 

controller, replacing the signal from the actuator model in each case (see 
figs. 5 and 6). This computation obviously cannot provide perfect information 
because of the inherent characteristics of the deadbands and actuator lags and 
will result in some coupling between the various controllers. Despite these 
imperfections and coupling, use of this technique met with considerable suc- 
cess as indicated in figure 15. It can be seen that the magnitude of the 
angular oscillations and the actuator motion are considerably reduced compared 
to the cases with no compensation or with model compensation. 

The integrated controller, using blade actuator signals fed back from the 
engine, was found to be tolerant of changes in actuator dynamics. An example 
is shown in figure 16, where it has been assumed that the time constant of the 
actuators will depend on whether the blade angle is increasing or decreasing. 
This difference in time constant is a possibility because of the presence of 
aerodynamic torques on the blades (see ref. 4). For the case shown in fig- 
ure 16, it is assumed that the time constant for increasing the blade angle 
will be 0.3 sec, three times the nominal value of 0.1 chosen for the study. 

The results show some increase in the angular oscillations, but the resulting 
pitch and roll oscillations due to the actuator activity remains well below 
1° amplitude. 

The effect of all changes to the blade actuator system studied are sum- 
marized in figure 17. These results are all for light turbulence and consist 
of maximum angular acceleration recorded during the time interval for the 
approach from 5 to 105 sec. The advantage of the feedback type mechanization 
in reducing the maximum pitch accelerations is clearly evident as is the 
decrease in performance when the time contant for increasing blade angles is 
tripled. A further doubling of the time constants for both increasing and 
decreasing blade angles roughly triples the maximum angular accelerations. 

With feedback compensation, changes in the deadband characteristics did 
not produce any large changes in the maximum angular acceleration. With no 
deadband compensation, a reduction in the deadband by a factor of 2 to 0.25° 
was found to reduce the motion by a factor of about 2. However, as indicated 
in figure 17, reducing the deadband caused only a small improvement for the 
compensated system. It was also found that different deadbands for each 
actuator, averaging to the nominal value of 0.5° (0.25°, 0.75°, and 0.50°), had 
negligible effect on maximum accelerations. 

In addition to the motion caused by actuator dynamics, the response of 
the compensated systems to increasing turbulence level was examined. These 
results are summarized in figure 18. If these results are extrapolated to 
zero turbulence, it can be seen that the increment contributed by the effects 
of fan blade pitch angle deadband are generally small compared to the overall 
effect of turbulence on angular accelerations. Further, the rate of increase 
of angular acceleration with increasing turbulence for no oeadband compensation 
or with feedback compensation differs little from that with zero deadband. 

This implies that the excellent response to turbulence of the state rate 
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feedback implicit model following controller (ref. 2) is retained in the pres- 
ence of fan actuator deadband. 

The flightpath performance showed no important variation for any of the 
controller configurations studied. For the lowest turbulence level 
(0.25 m/sec rms), the maximum attitude error was less than 0.11 m; for a tur- 
bulence level of 2 m/sec, this error was always less than 0.6 m. 


Effects of Temperature Rise Due to Engine Exhaust Reingestion 

The effects of reingestion were studied in the final stage of the approach 
to a hover. The commanded hover altitude was set at 4 m to ensure that the 
maximum effect of the temperature rise would be experienced. In practice it 
is unlikely that the approach trajectory would terminate at such a low alti- 
tude; instead hover would be established at a considerably higher altitude 
followed by a vertical descent to a landing. 

The temperature increase due to engine exhaust reingestion was estimated 
from the test results given in reference 7. The temperature rise index 
employed is shown in the upper part of figure 19. Data were not available 
from which to establish the temperature rise at a location corresponding to 
the inlet of the turboshaft engine mounted in the fuselage of the simulated 
aircraft. It was assumed that this temperature was identical to that at the 
inlet to the nose fan. This assumption was thought to be as realistic as any 
other and had the advantage of requiring a minimum of modification to the 
simulation. It should be noted that the data on which the temperature 
increases are based was obtained during a static test. Thus, no dynamic 
effects are included (e.g., horizontal or vertical velocities that exist for 
the trajectories studied). The results shown, therefore, represent a worst 
case situation since motion of the aircraft or any relative Ind would preclude 
the temperature from reaching steady state levels. 

Ground proximity also usually has an aerodynamic effect on propulsive 
lift. For this simulation, no such ground effect was assumed. 

The temperature variation at the inlet for a typical approach is shown in 
the lower part of figure 19. The curves have a somewhat different shape from 
those in the upper part of the figure since ingestion causes a positive feed- 
back whereby the temperature of the exhaust gases continues to Increase as the 
height decreases to about 4 m, causing a greater rate of Increase of inlet 
temperature. 

Simulation results were obtained for two approach types: the 3° curved 

approach and the 3° straight approach. In each case, the number two engine 
was assumed to have failed. The curved approach required a 0.61 m/sec 2 verti- 
cal acceleration to be reached during the last 8 sec to achieve hover. In the 
presence of an engine failure, this could not be achieved for gross weights 
greater than about 11,000 kg. The 3° straight approach does not require a 
deceleration normal to the flightpath, allowing a somewhat greater gross 
weight of 11,700 kg. The variation of altitude, net thrust, and turbine inlet 
temperature for each of the two approaches, with and without relngestlon, is 
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shown in figure 20. For these results wir.d and turbulence are assumed to be 
zero. As can be seen, for either approach, the reingestion temperature rise 
causes a decrease in the net thrust and a consequent loss in altitude of 2 to 
3 m before the controller could arrest the rate of sink. Also the turbine 
inlet temperature rise on the order of 120 K is imposed on the already elevated 
temperature as a consequence of the engine failure. The maximum temperature 
reaches about 1700 K or about 80 K in excess of the temperature stated in 
reference 6 as the 1 hr contingency rating. 


CONCLUDING REMARKS 


A nonreal-time digital simulation of a lift-fan VTOL aircraft with a 
detailed representation of the propulsion system has been employed to study 
interactions between the propulsion and flight control systems. In addition, 
the detailed propulsion system representation provided information useful for 
propulsion system design (e.g., temperature cycling data which may affect 
engine fatigue life). 

The results of the study indicated no deleterious effects of engine dynam- 
ics on the flight dynamics that could not be established using a much simpler 
engine model developed for real-time piloted simulation. The only deleterious 
effect observed on flight dynamics was caused by the deadband in the fan blade 
pitch actuators — an effect that can be included easily in the simpler propul- 
sion system representation. A technique was developed for alleviating the 
effects of this deadband through feedback of the actuator signals to the flight 
control system. 

Data were obtained documenting the response of critical internal engine 
parameters during an approach under autopilot control. Included are results 
that document the response to turbulence and engine failure and the effects of 
changes in ambient temperature and gross weight. These data are available on 
magnetic tape to agencies of the U.S. Government and their contractors. 
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TABLE 1.- AIRCRAFT DIMENSIONAL 
AND MASS DATA 


Wing span 

13.54 

m 


Wing area 

31.79 

m 2 


Empty weight 

10,886 

kg 


Fuel load 

2,041 

kg 


Moments of inertia, full 




fuel load 




lx 

30,973 

kg 

m 

l y 

53,971 

kg 

m 

Iz 

67,841 

kg 

m 

^xz 

3,022 

kg 

m 

Engine thrust 

151,057 

N 


Thrust -to-weight ratio — 




all engines operating 

1.19 
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TABLE 2.- LIST OF SIMULATED APPROACHES 


Approach 

number 

Turbulence 
m/sec, rms 

Engine 

failure 

number 

Fan actua- 
tor dead- 
band, deg 

Deadband 

compen- 

sation 

Other 

i 

1.00 

2 

c 

) 

— 


Initial altitude 









error - 20 m; initial 
localizer error = 20 m 

2 

.25 

No 



— 



3 

1.00 





— 



4 

2.00 




f 

— 



5 

.25 



0.5 

None 



6 

1.00 





None 



7 

2.00 





None 



8 

.25 





Model 


9 

1.00 





Model 


10 

2.00 





Model 


11 

.25 





Feedback 


12 

1.00 








13 

2.00 


f 






14 

.25 

2 






15 

1.00 

2 






16 

2.00 

2 






17 

2.00 


3 






18 

2.00 

2 





Ambient temperature for 









standard day, 15 °C 

19 

.25 

2 





Initial weight, 11,138 kg 

20 

2.00 

2 





Initial weight, 11,138 kg 

21 

.25 

No 





Fan blade actuator time 









constant for increasing 
blade angle, 0.3 sec 

22 

2.00 







Fan blade actuator time 









constant for increasing 
blade angle, 0.3 sec 

23 

.25 







Fan blade actuator time 









constant for increasing 
blade angle, 0.6 sec; 
decreasing 0.2 sec 

24 

2.00 







Fan blade actuator time 





1 




constant for Increasing 
blade angle, 0.6 sec; 





W 



decreasing 0.2 sec 

25 

.25 




.25 




26 

2.00 




.25 




27 

.25 



Variable 



Deadband #1, 0.25; 









#2, 0.75; #3, 0.50 

28 

2.00 



Variable 



Deadband #1, 0.25; 









#2, 0.75; #3, 0.50 
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TABLE 2.- CONCLUDED 


Note: The following conditions existed for all approaches, unless specified 

differently in the column designated "other." 


Flightpath 

Wind 

Initial gross weight 
Ambient temperature 
Fan blade actuator time 
constant 

Initial altitude and local- 
izer errors 

Altitude for engine failure 


Curved decelerating approach, tp = 90° — 
see figure 10 
20 knots, from ip - 60° 

12,960 kg 

32. 2°C (hot day) 

0. 1 sec 


0 

305 m 


TABLE 3.- LIST OF APPROACHES TO HOVER 


Approach 

number 

Type 

Gross 

weight 

Engine 

failure 

number 

Reingestion 

1 

3°-curved 

12,800 

No 

Yes 

2 



12,800 

No 

No 

3 



11,010 

4 

) 

No 

4 



11,010 



Yes 

5 



12,860 



No 

6 

3°-straight 

12,860 



Yes 

7 

3°-straight 

11,690 



Yes 

8 

3°-straight 

11,690 



No 


Note: All approaches are with blade actuator dead- 

band of 0.50 and with feedback compensation. No 
wind or turbulence. 
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Figure 1.- Sketch of simulated lift-fan aircraft. 
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Figure 2,- Block diagram of engine model 










Figure 3.- Attitude control schematic. 
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FIGURES 7 & 8 
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Figure 4.- Flight path control scheaetlc. 
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(T) HEIGHT CONTROL INPUT SCHEDULE 
(I) BLADE ANGLE BIAS SCHEDULE 

Figure 7.- Power lever system 
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Figure 8.- Fuel and engine control details. 
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Figure 9.- Fan blade pitch actuation details. 
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(a) Flightpath. 

Figure 10.- Airframe and engine data for a typical approach. Conditions: 
maximum gross weight, 20-knot wind 30° from flightpath headir.g (ip ■ 90®), 

1 m/sec rms turbulence, hot day, failure of no. 2 engine at 300 m altitude, 
no fan actuator deadband. 
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(b) Airframe variables. 
Figure 10.- Continued. 
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«8ure 10.- Continued. 
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(a) Pitch and roll attitude and acceleration. 


Figure 11.- Limit 
Conditions: 


cycle behavior induced by fan blade actuator deadbands. 
0.5° actuator deadband, 0.25 m/sec rms turbulence. 


31 


BLADE ANGLE COMMAND AND RESPONSE - FAN 3, deq 



(b) Details of typical fan blade angle command and response 
Figure 11.— Concluded. 
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Figure 12.- Effect of a 0.5 actuator de< 

freedota s; 
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Figure 13.- Deadband compensation technique for simple single degree-of- 

freedom system. 
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<i) Pitch and roll attitude and acceleration. 


Figure 14.- Lialt cycle behavior with model-type compensation for actuator 
deadband. Conditions: 0.5* actuator deadband, 0.25 m/sec ras turbulence. 
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TIME, sec 

(a) Pitch and roll attitude and acceleration. 

Figure 15.- Limit cycle behavior with feedback-type compensation for actuator 
deadbands. Conditions: 0.5° actuator deadband, 0.25 m/sec rms turbulence. 
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TIME, sec 


(a) Pitch and roll acceleration. 

Figure 16.- Response with actuator time constant for increasing blade angle 

tripled. 
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r/KN 3 , 




MAX. ANGULAR ACCELERATION, rad/*ec 2 


Figure 17.- Effect of deadband, compensation technique, and actuator time 
constants on maximum angular acceleration in time interval from 5 to 
105 sec. RMS turbulence * 0.25 m/sec. 
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Figure 18.- Effect of turbulence level on rms angular accelerations in time 

interval from 5 to 105 sec. 
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INLET TEMP. RISE INDEX 



ALTITUDE, m 


Figure 19.- Inlet temperature rise index and the calculated inlet temperature 

for a curved decelerating approach. 
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